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Abstract 

A magnetic origin for the honeycomb reconstruction of metallic surfaces with three-fold 
symmetry like Pb/Ge (111) is proposed. Assuming that the groundstate is an antiferro- 
magnet insulator over the triangular lattice of adatom sites (Pb) , we demonstrate that the 
former is simultaneously unstable to canting and to a structural distortion if the surface 
is soft enough. We therefore predict a net magnetization over the reconstructed surface at 
sufficiently low temperature. 
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The (111) surface of Ge covered with | of a monolayer of Pb (or Sn) adatoms displays a 
striking reconstruction upon cooling. 1-3 At room temperature the surface is metallic, with 
the adatoms forming a triangular lattice. However, at low temperatures, T < T* = 250 K, 
a honeycomb relief forms over the triangular adatom structure. This transition from a 
x \f?> surface structure at ambient temperature to a reconstructed 3x3 surface at 
low temperature is observed directly via scanning electron microscopy (STM) and by low- 
energy electron diffraction (LEED) techniques. The reconstruction is also reflected in 
the electronic structure of the surface. In particular, both angle-resolved photo emission 
spectra (ARPES) and electron energy loss spectra (EELS) find that a pseudogap opens 
at the Fermi surface of the surface conduction band at temperatures below T*. 4 ' 5 Such 
measurements suggest that the reconstructed surface has an insulating groundstate. 

It was originally proposed that the 3x3 surface reconstruction is due to a charge- 
density wave (CDW) instability within the metallic \/3 x \/3 surface. 1 The absence of good 
nesting properties of the Fermi surface at room temperature 2-5 indicates, however, that 
this is unlikely. To account for the phenomenon that remains, we propose instead that the 
reconstructed surface is an antiferromagnetic Mott insulator 6 ' composed of localized spin- 
1/2 moments associated with the dangling chemical bond left at each adatom site. Both 
the narrowness of the surface conduction band 4 ' 8 and the presence of the low-temperature 
pseudo-gap feature 1 ' 5 make this idea plausible. In the limit of large on-site repulsion, we 
are then left with an antiferromagnetic Heisenberg model description for the electronic 
spins that remain over the sites of the triangular lattice of adatoms. A classical analysis 
of this model in combination with magnetostriction phenomenology 9 leads us to conclude 
that the magnetic frustration intrinsic to antiferromagnetism on a triangular lattice drives 
a soft enough v^3 x v^3 surface to buckle or dilate with the periodicity of a honeycomb relief. 
We therefore predict that the low-temperature 3x3 surface reconstruction has a struc- 
tural component. The original LEED evidence for the reconstruction of the Pb/Ge (111) 
surface 1 as well as recent surface X-ray diffraction studies 10 ' 11 support this point of view. 
Coincident with the honeycomb distortion is a canted antiferromagnetic spin arrangement 
with reduced frustration (see Fig. 1). We therefore also predict a net magnetization over 
the reconstructed surface at low enough temperatures. 

We begin the theoretical discussion at zero temperature, where a Mott insulator 
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state 6 ' 7 is assumed for the reconstructed metallic surface. 1 Specifically, suppose that a 

— * 

localized spin s = ^ moment, Si, lies at each site, i, of the triangular lattice coincident 
with the Pb (or Sn) adatoms. The former represents the spin associated with the dangling 
chemical bond on the Pb/Ge (111) surface. The simplest description of this system is an 
antiferromagnetic Heisenberg model 

H = ^(JijSi ■ S, + 6J i:j S?S*) -h-J2$i (!) 
(y) * 

summed over nearest-neighbor bonds (ij) on the triangular lattice. Here, J t j > denotes 

the exchange coupling constant over such bonds, while h denotes the external magnetic 

field. Also, the parameter 6 is a measure of the magnetic anisotropy along the z-axis 

perpendicular to the surface that originates from the spin-orbit interaction. We now take 

the key step in the theory by supposing that a honeycomb lattice distortion accompanies 

the Mott insulator state (1), and that this results in two different exchange coupling 

constants, J and J', over the undistorted and distorted bonds, respectively (see Fig. 1). 

The buckled (dilated) lattice is achieved by raising (expanding) or lowering (contracting) 

the regions in the vicinity of the honeycomb (a and b sites) with respect to those regions 

in the vicinity of the centers (c sites) of the honeycomb. The phenomenology 

Jij = Jo + Ji(Vi + Vj) ( 2 ) 

for the exchange coupling constant of the distorted lattice is compatible with these assign- 
ments. Here J\ is the first order magnetostriction coefficient 9 with respect to the scalar 
field rji, which can represent either perpendicular displacements, u z , or lateral dilations, 
d x Ux + dyUy, of the (triangular lattice) surface in the vicinity of adatom i. Equations (1) 
and (2) represent the essential components of the theory. 

We now analyze the isotropic (5 = 0) antiferromagnetic model (1) in the classical 
limit, s — > oo, in which case the spin operator Si reduces to a 3-vector of magnitude s. 
The ground state for the honeycomb arrangement of exchange coupling constants J a b = J 
and J ac = J' = Jbc indicated in Fig. 1 is a canted three-fold symmetric spin configuration 
over the three sublattices a, 6, and c, of the honeycomb-distorted lattice. This canted 
antiferromagnet has a net magnetization 

S = l -{S a + S b + S c ) = 1 (1 - 2 cos 9)S C (3) 
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aligned parallel to the sublattice magnetization, S c , of the c sites, where 9 represents the 
half angle in between spins S a and Sb (see Fig. 1). Since this macroscopic moment must 
be aligned along the external magnetic field, h, we obtain an energy per site 

E M /N = s 2 Jcos(29) + 2s 2 J'cos(tt -9)- E Z \S\ (4) 

for the canted antiferromagnet, where E z oc \h\ denotes the Zeeman energy splitting. If 
we define A J = J' — J, then the minimum energy occurs at a half angle 9q such that 
cos# = [s 2 J' + l(sgnAJ)sE z ]/2s 2 J, with energy 

9 [sJ' + \{sgnAJ)E z ] 2 1 , a \ 
E M /N = -s 2 J- 1 - '-j£L + -(sgn AJ)sE z , (5) 



and with net magnetization 



1/AJ sgnAJ E z \^ 



)S C . (6) 



31 J 3 sJ 

Notice that both an external surface strain [see Eq. (2) and ref. 9] and an external 
magnetic field increases canting, and the saturated magnetization, S, as a result. 

Above, we have shown how a honeycomb distortion of the frustrated antiferromagnet 
on the triangular lattice lowers the magnetic energy via magnetostriction. But is such a 
canted antiferromagnetic state energetically favorable with respect to the distorted lattice 
itself? We return to the magnetostriction phenomenology (2) to address this question. 
Suppose that the triangular lattice has a honeycomb distortion, with lateral dilations 
or perpendicular displacements r] a = rji, at the honeycomb, and with lateral dilations or 
perpendicular displacements rj c at the centers. Macroscopic stability of the surface requires 
that 

Va+Vb + Vc = 0. (7) 

This fact coupled with expression (5) for the magnetic energy in the absence of magnetic 
field, in addition to the magnetostriction phenomenology (2), yields the new expression 

E M /N = -s 2 [J Jo + J 2 (Va - Vc) 2 ] (8) 
that is valid up to quadratic order in the displacements, where 



is the effective second order magnetostriction coefficient of the distorted lattice. Therefore, 
if the bare energy for a lattice distortion has the effective form E L = \k Q Yl^iVi ~~ Vj) 2 i 
then the system as a whole, E = El + Em, is unstable to a honeycomb- type buckling or 
dilation for soft surfaces such that 

s 2 J 2 > 2k . (10) 

Here, ko denotes the relevant effective spring constant for the case of a buckling distortion. 
We therefore observe that the frustrated antiferromagnetism that exists on the triangular 
lattice of spins drives the honeycomb distortion in order to improve magnetic energy on 
the soft surface (10). Last, Eqs. (4)-(6) indicate that this effect is only enhancedbj the 
application of an external magnetic field! 

Yet what effect do the quantum fluctuations that are connected with the localized spin- 
1/2 moments have on this honeycomb- type lattice instability? It is known that frustration 
generally decreases the spin stiffness of a spin-1/2 antiferromagnet. 12 In our notation, the 
spin stiffness is essentially the prefactor s 2 in expression (8) for the magnetic energy. This 
means that quantum corrections renormalize down s 2 to a function s 2 R of frustration that 
has a minimum at J = J', in which case the frustration is at a maximum. Hence, quantum 
corrections can only enhance the honeycomb lattice instability (10) driven by frustration 
in the antiferromagnet (relative to the undistorted triangular lattice [see Eq. (8)]). 

We have seen how a canted antiferromagnetic spin arrangement over the triangular 
lattice can induce a buckling or dilation instability with a periodicity of the honeycomb 
type. Strictly speaking, the effect requires long range order among the magnetic moments. 
A "large- N n analysis of the continuum energy functional 

E M = \ Ps j d 2 r[(Vm) 2 ± ^ 2 m 2 z ] (11) 

with respect to the normalized sublattice magnetization (|m| = 1) of the effective classical 
Heisenberg model (1) indicates that this is indeed the case at temperatures below a critical 
temperature, 13 ksT c = 27r / o s /ln|5| _1 . Here, p s = s 2 R J denotes the relevant spin stiffness, 13 
while 5' = a' 2 /£ 2 is a measure of the magnetic anisotropy along the z-axis perpendicular 
to the 3x3 surface with lattice constant a'. The anisotropy 5' is assumed to be small. 
A net macroscopic moment is then predicted to exist over the (3 x 3) surface at such 
low temperatures, T < T c . At high temperatures, T > T c , the magnetic correlation 
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length, £ M {T), is finite. Clearly, the honeycomb lattice instability (10) cannot occur if the 
magnetic correlation length is short compared to the lattice constant. This suggests the 
(implicit) definition £m(T*) = a' for the paramagnetic cross-over temperature, /c^T* ~ 
2np s . At intermediate temperatures T c < T < T*, we expect to have fluctuating domains 
of ferrimagnetic order, each with dimensions of order £m(T). A short-range honeycomb 
lattice instability (10) is then possible. In such case, we expect a diffuse Bragg diffraction 
pattern corresponding to the 3x3 surface. 

In the regime of high temperatures, T > T* ~ J/ks, our simple model [Eqs. (1) and 
(2)] thus predicts a paramagnetic insulator state over an undistorted triangular lattice. 
This prediction, clashes, however, with experimental evidence for the existence of a metallic 
surface at ambient temperatures above the reconstruction temperature, T c . The electronic 
charge fluctuations that are suppressed in the Heisenberg model (1) must therefore be 
included. We resort to the nearest-neighbor Hubbard model 

H = - *« c t c jo- + u i n ^ n n ( 12 ) 

on the triangular lattice at half filling to address this problem. Here, c ia denotes the 
annihilation operator for an electron of spin a at site i, while is the corresponding 
occupation number. In the limit of strong on-site repulsion, Ui — > oo, we recover the 
antiferromagnetic insulator state (1) already discussed, with exchange coupling constants 6 

J ij = 2t 2 ij (U~ 1 + U- 1 ). (13) 

Let us consider the weak-coupling limit Ui = instead, and suppose that a honeycomb lat- 
tice distortion is "frozen in" . The electron-phonon interaction then results in two different 
nearest-neighbor hopping amplitudes, t a b = t and t ac = t' = tb c , in between the sublattices 
(see Fig. 1). Assuming Bloch waves of momentum k for the one-electron states over each 
sublattice, we obtain a corresponding spectrum that satisfies the characteristic equation 

= £ |- (2t' 2 + t 2 )\a % \\ + t'H{a\ + af\ (14) 

where the amplitude 

a- = e itdl + e lts2 + e^ 3 (15) 
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is a sum over the (overcomplete) basis, a\ = ax and a 2) 3 = (a/2)(— x ± 3 1 / 2 ?/), of the 
triangular lattice (see Fig. 1). This system has two noteworthy limits: (i) On the undis- 
torted triangular lattice, t = t' , we have a single band of surface conduction electrons that 
disperses as 

e £ = + «£)• ( 16 ) 

The density of states per site at the Fermi surface for the half-filled band is N(0) ~ W~ l , 
where W = 9t gives the bandwidth. The Fermi surface itself at half filling is nearly 
circular and evidently shows no nesting (see Fig. 2, at T > T c ). (ii) In the case of an 
extreme honeycomb distortion, t <C t' or t' <C t, we have two graphene-type bands that 
are separated by a small gap at the vertices (if 3x3 ) of the Brillouin zone boundaries that 
correspond to the honeycomb. Within this gap lies a narrow band, 

e s = f(a| + aj 3 )/|a £ | 2 , (17) 

that crosses the Fermi surface at half filling, with an effective hopping amplitude i = \t 
and i = t' 2 /t in the respective limits t — > and t' — > 0. This band has a global maximum 
at the center of the Brillouin zone (r3 X 3) and a global minimum at the points M3 X 3 that 
lie in between neighboring -K3X3 points. The Fermi surface at half filling for this low- 
temperature (T < T c ) band structure is displayed in Fig. 2. The combined effect of a 
narrow bandwidth, W = 8t, and of the "corners" present in the Fermi surface result in a 
relatively large density of states, N(0) ~ VF -1 , there. 

In light of these observations, we propose the following mechanism that accounts for 
the metallicity of the unconstructed surface. At high temperatures T > T c , the electronic 
states at the surface are described by a half-filled nearest-neighbor Hubbard model (12) over 
the triangular lattice, t = with an on-site repulsion U below the critical one, U c ~ W, at 
the Mott transition. 7 This is consistent with a metallic state. At low temperatures, T < T c , 
the surface buckles or dilates (t 7^ t') in order to lower its magnetic energy (J 7^ J'). 
In particular, the honeycomb distortion of the triangular lattice of adatoms results in a 
narrow band (17) that crosses the Fermi surface through the effects of the electron-phonon 
interaction. The presence of an insulating magnetic groundstate at zero temperature then 
requires that the on-site repulsion U now be greater than the new critical one, U c ~ W, at 
the Mott transition. 7 This means, in rough terms, that the on-site repulsive energy must 
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lie inside the window W < U < W of parameter space. Last, the Fermi surface of the 
reconstructed metallic surface (see Fig. 2, at T < T c ) should be unstable to the formation 
of an SDW due to the relatively large density of states connected with the narrow low- 
temperature band (17). 14 Since the electron-phonon interaction induces magnetostriction 
(2) through the identity (13), the previous strong-coupling analysis indicates that this 
SDW state is precisely the canted antiferromagnetic arrangement of spins depicted in Fig. 
1. We then notably predict that the reconstructed surface has a net magnetization (6) at 
low temperatures, T < T c . An analysis of the Hubbard model (12) within the mean-field 
approximation indicates that such a net moment causes the conduction bands (17) of the 
reconstructed surface to split in a manner similar to Fe. 8 

Yet how does the above proposal compare with experiment? First, the presence of an 
SDW gap at the Fermi surface is consistent with evidence that points towards the existence 
of an insulating state on the reconstructed Pb/Ge (111) surface at low temperature: e.g., 
STM images, EELS and surface photo-emission spectra. 1 ' 4 Second, the band structure 
characteristic of a honeycomb distortion of the triangular lattice (17) is consistent with 
recent ARPES studies of the reconstructed Pb/Ge (111) surface 5 that see the gap expected 
by zone folding at the M^xs points of the Brillouin zone. Even more striking is the fact 
that such ARPES experiments observe iron-type splitting of surface conduction electrons 
bands! This of course is consistent with the presence of a net magnetic moment, 8 which 
we propose is due to the canted antiferromagnetic Mott insulator state shown in Fig. 1. 
Last, both LEED and recent X-ray diffraction studies find evidence for a periodic (3 x 3) 
structural distortion of the reconstructed Pb/Ge (111) surface that is consistent with the 
honeycomb distortion proposed here. 1 ' 10 ' 11 

Contrary to expectations based on the width of the surface conduction electron band, 2 
however, the 3x3 reconstruction presented by the Sn/Ge (111) surface appears at temper- 
atures below T* = 210 K, which is appreciably lower than that corresponding to the Pb/Ge 
(111) surface. And unlike the latter, the Sn/Ge (111) surface appears to remains metal- 
lic down to a temperature T = 110 K. 2 ' 3 We speculate that such discrepancies between 
the Pb and Sn adatoms surfaces originate from the elemental differences in the spin-orbit 
coupling. The relatively weak spin-orbit coupling present in the Sn adatom results in a 
smaller magnetic anisotropy 5 [see Eqs. (1) and (11)], and thus in both a lower magnetic 
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rigidity, /csT*, and critical temperature, T c . 
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Figure Caption 



Fig. 1 Shown is the canted antiferromagnetic spin arrangement over the 3x3 lattice. 

Fig. 2 The Fermi surface (FS) for both the low-temperature (17) and the high-temperature 

(16) band structures at half filling are displayed. Also displayed are the corresponding 

Brillouin zone (BZ) boundaries. 
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